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Synergy Dynamics of Vortices and Solitons in an Atomic Bose-Einstein Condensate 

Excited by an Oscillating Potential 
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The hydrodynamics of quantized vortices and solitons in an atomic Bose-Einstein condensate 
excited by an oscillating potential are studied by numerically solving the two-dimensional Gross- 
Pitaevskii equation. The osciUating potential keeps nucleating vortex dipoles, whose impulses al- 
ternatively change their direction synchronously with the oscillation of the potential. This leads 
to synergy dynamics of vortices and solitons which have not been previously reported in quantum 
fluids. 

PACS numbers: 67.85.De,03.75.Lm,67.25.dk,47.37.-|-q 



I. INTRODUCTION 

Topological defects are key concepts in general physics 
[1]. Quantum condensed systems are ideal for studying 
topological defects. Quantized vortices and solitons have 
been thoroughly characterized in superfluid ^He [2] and 
■^He [3, 4]. Modern research developments on quantum 
turbulence (QT) is still based on understanding dynam- 
ics of quantized vortices [5]. An atomic Bose-Einstein 
condensate (BEC) can be used to investigate quantum 
hydrodynamics. Atomic BECs have several advantages 
over superfluid helium. In particular, modern optical 
techniques enable control of the condensate and direct 
visualization of topological defects such as vortices and 
solitons. Actually many important works have been per- 
formed on quantized vortices these years [6]. Recently 
Henn et al. made and observed QT in a BEC by intro- 
ducing an external oscillatory perturbation of the trap- 
ping potential [7]. 

Here we numerically address the response of a BEC to 
an oscillating repulsive potential. Vibrating structures 
such as spheres, grids, and wires are used in superfluid 
''He and ^He to create QT [5, 8]. Despite the differences 
between these structures, the experiments show surpris- 
ingly similar behavior. We apply this methodology for 
atomic BECs for the first time. 

A few works on oscillating potentials in atomic BECs 
have already been numerically and experimentally per- 
formed [9-11]. Dissipation which works above some crit- 
ical velocity was investigated theoretically by Jackson et 
al. [9] and experimentally by Raman et al. [10]. How- 
ever, they did not observe the dynamics of vortices and 
solitons reported in this work. Recently, Neely et al. ob- 
served the formation and dynamics of vortex dipoles by 
forcing superflow around a repulsive Gaussian potential 
within an oblate BEC [12]. 

Studying quantum hydrodynamics in an atomic BEC 
subject to an oscillating potential should open up a re- 
search area different from helium and other BEC cases 
for the following reasons. First, since an atomic BEC is 
a clean system free of remnant defects and impurities, 
it is possible to study the intrinsic nucleation of topo- 



logical defects. Second, the oscillating potential leads to 
synergy dynamics of quantized vortices and solitons that 
have not been previously observed. Linear motion of the 
potential results in vortex dipoles with fixed charge [12], 
whereas oscillatory motion leads to vortex dipoles with 
alternating charges, causing rearrangement of dipoles and 
a metamorphosis between vortices and solitons. Third, 
the oscillation introduces another important parameter, 
namely frequency, into quantum hydrodynamics. Nucle- 
ation of defects and the transition to QT depend on fre- 
quency, which should be investigated in future. Fourth, 
this work develops a powerful new method for making 
QT in a trapped BEC in addition to other known ones 
[7, 13, 14]. Eventually the dynamics of vortices and soli- 
tons can be visualized in atomic BECs, enabling a direct 
comparison between experiments and theoretical or nu- 
merical results. 

In this work we investigate the quantum hydrodynam- 
ics of a trapped BEC with an oscillating potential by 
numerically solving the two-dimensional Gross-Pitaevskii 
(GP) equation. In Section II the model and numerical 
calculation are described. Section HI presents the whole 
dynamics of vortices obtained by our numerical calcula- 
tion. Next we explain elementary processes related to 
vortices and solitons in Section IV. The oscillating po- 
tential may heat the condensate, which is estimated in 
Section V. Finally we summarize our work in Section VI. 



II. THE MODEL AND NUMERICAL 
CALCULATION 

We consider a dilute atomic BEC, assuming that the 
condensate is pancake shaped. This system is well de- 
scribed by a macroscopic wavefunction ip obeying the GP 
equation 
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where m is the particle mass, V is the potential, and g 
is an interaction parameter for the two-dimensional case. 
The wavefunction ip is normalized by the total particle 
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FIG. 1: (color online) Migration of vortices. Density profile at (a)t=Os, (b)f=0.0345s, (c) t=0.0591s, (d) t=0.0760s, (e) 
i=0. 0845s, (f)i=0.f49s, (g)i=0.21f s, and (h)f=0.240s. The x and y dimensions of the images are 140 ^m and 32.5^m. The 
symbol + (-) denotes a vortex with clockwise (counter-clockwise) circulation. There are two kinds of holes. The central large 
hole is the Gaussian potential, while the other small holes are vortices. 



number N . Suppose the condensate is confined by a har- 
monic potential 14 and penetrated by a Gaussian poten- 
tial Vg, so that V = Vh+Vb where 14 = im(a;^a;2-hw^y^) 
and Vg = Vbexp[— ((x — xolt))"^ + y'^)/(P]. Here xo(i) is 
the x-coordinate of the center of the Gaussian poten- 
tial and d is its radius. Oscillate the Gaussian potential 
as a;o(i) = esin(ajt). We use a dimensionless form of 
Eq. (1) for numerical calculations. Space and time are 
normalized by h/ y/2mgnQ and h/gno, where no is the 
density near the center of the condensate. Choose pa- 
rameters g = 4.19 X 10"''^ J/m^, m = 1.42 x 10"^^ kg, 
iV = 6.6 X 10**, u}^ — 2n X 5/s, u}y == 2it x 25 /s, 
d = 0.6 /xm, e ~ 7 fim, u = 100 /s and Vq — QQgriQ. 
We use Crank-Nicholson method to perform numerical 
calculations without dissipation and noise. For the simu- 
lation, space in the x and y directions is discretized into 
2048x640 bins. 



III. WHOLE DYNAMICS 



ternately change direction. They reconnect with each 
other to make new vortex pairs, leaving the potential 
in Fig. 1(c). This phenomenon is not observed for the 
case of uniform motion of the potential, but only for an 
oscillating potential. Reaching the surface in Fig. 1(d), 
the vortex pairs interact with ghost vortices, which are 
vortices in the low-density region and described later in 
IV. A. Then the vortices head toward the bow of the con- 
densate along the surface in Fig. 1(e) and (f). A vortex 
coming up from the left side reaches the bow to meet 
one from the right side in Fig. 1(g), thus making a new 
vortex pair. Finally, the pair comes back to the center of 
the condensate in Fig. 1(h). Thus the vortices nucleated 
by the potential enjoy a lengthly migration in the "sea" 
of BEC; the vortices are nucleated from the potential, 
reconnect, move away from it, reach the surface, head 
toward the bow and come back to the center. In the fol- 
lowing, we illustrate elementary processes related to the 
synergy dynamics. 



The dynamics of vortices in which they experience a 
lengthy migration are shown in Fig. 1. Following the des- 
tiny of vortices nucleated by the oscillating potential en- 
ables us to survey their dynamics. The initial state in 
the static Gaussian potential in Fig. 1(a) is obtained by 
an imaginary time step of the GP equation. A vortex 
pair is nucleated behind the Gaussian potential in Fig. 
1(b) as the potential starts to move. Then the oscillat- 
ing potential nucleates vortex pairs whose impulses al- 



IV. ELEMENTARY PROCESSES 



In this section we discuss elementary processes, which 
are nucleation of vortices, reconnection of vortex pairs 
near the potential, divorce of vortex pair at the surface, 
nucleation of solitons and collapse of solitons. These pro- 
cesses occur in the dynamics of Fig. 1. 




FIG. 2: (color online) Vortex nucleation by the potential: (a) 
i=0.00253s, (b) t=0.00490s, and (c) t=0.0118s. The upper 
and lower panels are density and phase profiles. The x and 
y dimensions of each image are 13 /im and 8/^m. The sym- 
bol + (-) denotes a vortex with clockwise (counter-clockwise) 
circulation. Black arrows in the density profile indicate the 
direction of motion of the potential, while the white arrows in 
the phase profile show the backfiow. The value of the phase 
varies from -tt (white) to tt (black). 
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FIG. 3: (color online) Reconnection of vortices near the po- 
tential: Density profile at (a)f=0.0507s, (b)f=0.0551s, (c) 
i=0. 0571s, (d) t=0. 0596 s. The x and y dimensions of each 
image are 39 fim and 24 /xm. The symbol + (-) denotes a 
vortex with clockwise (counter-clockwise) circulation. Black 
arrows indicate the direction of motion of the potential. The 
dashed arrows indicate the direction of the motion for vortices 
in (d). 



A. Nucleation of vortices 

A key issue is to understand how quantized vortices are 
nucleated. For a condensate to have vortices, it should 
bring some seeds of the topological defects from some- 
where. This situation is reminiscent of the formation of 
a vortex lattice in a rotating BEC [15]. Rotation nu- 
cleates "ghost vortices" in the low-density region at the 
outskirts of the condensate, and they are dragged into the 
interior to become the usual vortices having the conden- 
sate density. Our condensate has to do a similar thing. 
Ghost vortices can appear in two low-density regions in 
our system; one is inside the oscillating potential and the 
other is at the outskirts of the condensate. As soon as 
the potential starts to oscillate, a pair of ghost vortices is 
nucleated inside the potential and the condensate surface 
is filled with them. 

Vortex nucleation induced by the potential is shown in 
Fig. 2. It follows that ghost vortices inside the potential 
cause vortex nucleation. We consider nucleation of both 
kinds of vortices in this work. There are two critical ve- 
locities; Vci for a pair of ghost vortices and Uc2 for a pair 
of vortices with Vd < Wc2- When the velocity Vp = eu oi 
the potential is smaller than t^ci, vortices are not nucle- 
ated in the condensate. When Vp exceeds Wci, a pair of 
ghost vortices is nucleated inside the Gaussian potential. 
If Up is smaller than i;c2, however, those ghost vortices are 
annihilated. The repeated nucleation and annihilation of 
pairs of ghost vortices are characteristic of the oscillating 
case with Vd < Vp < Wc2. The observable critical velocity 
is Uc2 ~ 580/im/s for the present case, but tic2 strongly 
depends on the frequency w of the oscillation, which is 
an important target in future. The critical velocity Vc2 is 
smaller than the sound speed 1.33 x 10^/im/s, which is 
reported by Jackson et al. [9] . When Vp exceeds z;c2 , vor- 
tices are nucleated. These processes are shown in Fig. 2. 
A velocity field with backfiow is induced by the potential 



in Fig. 2(a). Then a pair of ghost vortices is nucleated in- 
side the potential in Fig. 2(b). Finally, the ghost vortices 
move away from it and a vortex pair appears in Fig. 2(c). 



B. Reconnection of vortex pairs near the potential 

Figure 1(c) shows a phenomenon characterized by the 
oscillating case after vortex nucleation, as detailed in Fig. 
3(a) to (d). After a vortex pair is nucleated behind the 
potential, it has an impulse in the direction of the veloc- 
ity of the potential. Hence the pair follows the potential 
in Fig. 3(a). Thereafter the oscillating potential changes 
the direction of its velocity in Fig. 3(b) and it makes a 
new ghost vortex pair whose direction of impulse is op- 
posite to that of the old pair. Then the ghost pair gets 
away from the potential, becoming a new vortex pair. 
The pair immediately reconnects with the old pair in Fig. 
3(c). This process produces new vortex pairs leaving the 
potential, and the pairs go toward the surface of the con- 
densate in Fig. 3(d). 



C. Divorce of vortex pairs at the surface 

The vortex pair escaping from the potential in Fig. 
3(d) reaches the surface and interacts with many fiuctu- 
ating ghost vortices there. The pair divorces, and the 
two resulting vortices move in opposite directions along 
the surface. The motion is caused by the following mech- 
anism. The normal component of the velocity field at the 
surface is suppressed when vortices are close to the sur- 
face. This makes the velocity field parallel to the surface, 
which moves the vortices along the surface. This motion 
may be understood shortly by using the idea of image 
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FIG. 4: (color online) Nucleation of solitons: (a) f=0.0929s, (b)f=0.106s, (c)t=0.111s, (d)t=0.114s. The panels (a)~(d) and 
(e)~(li) show the density(upper) and phase(lower) profile. The x and y dimensions of (a)~(f) and (g), (h) are 39 /im x 24 /im 
and 9 fim x 6 fim. The phase profiles for (c) and (d) are respectively (e) and (f). The figures (g) and (h) are the enlarged 
figures corresponding to the square boxes in (e) and (f). The symbol + (-) denotes a vortex with clockwise (counter-clockwise) 
circulation. Black arrows indicate the direction of motion of the potential in (a)~(d), but it shows the velocity field in (h). In 
the process from (a) to (b), reconnection of vortices occurs, then two new vortex pairs, which will become solitons, are created. 
These pairs are depicted by closed loops with dashed lines. The closed loops with solid lines show the solitons. The dashed 
arrows indicate the direction of the motion for vortices in (c), and one for solitons in (d). The value of the phase varies from 
-TT (white) to n (black). 



vortices [16], though the idea assumes that the boundary 
is a soHd waU and this is not exactly the case. 







FIG. 5: (color online) Collapse of solitons: Density profile at 
(a)i=0.126s, (b)i=0.131s, (c)i=0.136s, (d)t=0.138s. The 
X and y dimensions of each image are 39 /im and 24 /.im. 
The symbol + (-) denotes a vortex with clockwise (counter- 
clockwise) circulation. Black arrows indicate the direction of 
motion of the potential. The closed loops with solid lines show 
the solitons. The symbols © and from (b) to (d) denote 
the vortices created through the collapse of the soliton. The 
dashed arrows indicate the direction of the motion for solitons 
in (a). 



D. Nucleation of solitons 



The most interesting phenomenon in our research is 
the synergy dynamics between solitons and vortices. Nu- 
cleation, collision, and collapse of solitons are observed. 
These occur after Fig. 1(e). A soliton can be identified 
by its low density and phase jump of about tt. Figure 
4(a) through (h) show the nucleation of a soliton, which 
is peculiar to the dynamics by an oscillating potential. 
After the process in Fig. 3(a) to (d) occurs twice in both 
the right and left regions, the configuration of Fig. 4(a) is 
obtained. The vortex pairs in the right region in Fig. 4(a) 
get near the surface, and reconnect with other vortices in 
Fig. 4(b). The vortex pairs are depicted by closed loops 
with dashed lines. The pairs head toward the center of 
the condensate in Fig. 4(c). Subsequently, annihilation of 
a pair occurs as it reaches the center, becoming a soliton 
in Fig. 4(d). The change from a vortex dipole to a soliton 
is related to the nucleation of rarefaction pulse [17, 18]. 

One may think that the edges of the solitons have 
quantized vortices so that the low density regions inside 
the closed loops in Fig. 4(d) are merely vortex dipoles. 
However, there is a clear distinction between a vortex 
dipole and a soliton, which is obviously seen in the phase 
profiles (Fig. 4(e), (h)) corresponding to Fig. 4(c) and (d). 
Figure 4(e) and (f) show the phase profiles before and af- 
ter the annihilation of the dipoles. Their enlarged phase 
profiles are Fig. 4(g) and (h), which show respectively the 
vortex dipole and the soliton. Figure 4(g) shows appar- 
ently that the dipole has the topological defects in the 
phase profile, while the soliton in Fig. 4(h) does not have 



such defects and the velocity field around the edges does 
not have quantized circulation as seen by the arrows re- 
ferring to superflow. There remains the velocity field 
related to quantized vortices outside the edges, whereas 
the velocity field which tends to cancel the outside flow 
are induced in the low density region between the edges. 
Thus we can identify the soliton by checking the density 
and phase profile. 



E. Collapse of solitons 

The collapse of the solitons nucleates vortex pairs in 
(a) to (d) of Fig. 5. The solitons in Fig. 4(d) have opposite 
impulse, so that they collide with each other. However, 
they pass through each other without changing shape, 
in accord with the nature of solitons. The solitons then 
move toward the surface in Fig. 5(a). A two-dimensional 
soliton is unstable, and collapse of the solitons occurs at 
the surface. As a result, a soliton decays into a vortex 
pair at the surface, as shown in Fig. 5(b). These vortex 
pairs are depicted by closed loops with solid lines. Imme- 
diately after the collapse, the vortex pairs interact with 
ghost vortices in Fig. 5(c). Then, the vortices move to- 
ward the bow of the condensate along the surface in Fig. 
5(d) because of the mechanism described in IV. C. This 
process repeats as long as the Gaussian potential oscil- 
lates. This sort of transformation between vortex dipoles 
and solitons has been reported by Huang et al. [19]. 



The change is found to be only 1 nK even by the time of 
the last image in Fig. 1. Hence the heating is negligible, 
and use of the GP model is valid. 



VI. CONCLUSION 

We have performed numerical calculations of the two- 
dimensional GP equation to investigate the dynamics of 
vortices and solitons in a trapped BEG induced by an os- 
cillating potential. This paper reveals the dynamics char- 
acterized by oscillation. It is essential that the oscillat- 
ing potential makes vortex dipoles with different charges, 
which results in the synergy dynamics between vortices 
and solitons. The parameters used in these calculations 
are appropriate to experiments on atomic BEGs. We 
confirm the dynamics not only for the single fixed geom- 
etry but also for similar sets of parameters and geome- 
try. Hence, the dynamics obtained in our study should 
be experimentally observable. However, the dynamics by 
the oscillating potential should be dependent on the fre- 
quency and amplitude of the oscillation. Preliminary cal- 
culations indicate that a potential with lower amplitude 
creates more solitons, in contrast to the present result. In 
addition, the aspect ratio of the condensate shape should 
affect the dynamics. The systematic studies of how the 
dynamics depends on the parameters or geometry are re- 
ported elsewhere. 



V. HEATING OF THE CONDENSATE 
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